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Abstract 

The recent observations of sizable transverse fractions oi B ^ (pK* may hint for the existence of 
new physics. We analyze all possible new-physics four-quark operators and find that two classes of 
new-physics operators could offer resolutions to the B (f)K* polarization anomaly. The operators 
in the first class have structures (1 — 75) (1 — 75),cr(l — 75) (8'(t(1 — 75), and in the second class 
(1 +75) 0(1 +75), cr(l +75)0 (7(1+75). For each class, the new physics effects can be lumped into a 
single parameter. Two possible experimental results of polarization phases, aig{A±) — arg(A||) ~ vr 
or 0, originating from the phase ambiguity in data, could be separately accounted for by our 
two new-physics scenarios: the first (second) scenario with the first (second) class new-physics 
operators. The consistency between the data and our new physics analysis, suggests a small new- 
physics weak phase, together with a large(r) strong phase. We obtain sizable transverse fractions 
^11 II + ~ ^00) in accordance with the observations. We find A|||| ~ 0.8A^_|_ in the first scenario 
but A|||| > A_L^ in the second scenario. We discuss the impact of the new-physics weak phase on 
observations. 
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I. INTRODUCTION 



The studies for two-body charmless B decays have raised a lot of interests among the particle 
physics community. Recently the BABAR and BELLE collaborations have presented important 
results for the B meson decaying to a pair of light vector mesons (with V = 4>,p, or K*) P, 
1^ 1^ Q, IB; 1^ • This immediately surges a considerable amount of theoretical attentions to study 
nonperturbative features or to look for the possibility of having new physics (NP) in order to explain 
several dis crep ancies between the data and the Standard Model (SM) based calculations QjflS) 



several aiscreps 

HQ 0,0. 



From the existing SM calculations for the charmless B — > VV modes, it is known that 
the amplitude, Hqq (~ 0{1)), with two vector mesons in the longitudinal polarization state 
is much greater than those in transverse polarization states, since the latter are found to be 
H— ~ 0{l/mb),H++ ~ 0{l/ml) (or A\\ ~ ~ 0{l/mb) in the transversity basis) Q]. For the 

B meson decays, the relation for different helicity amplitudes is modified as Hqq : : H ~ 

0(1) : 0(l/?7ifc) : 0{l/m^). Nevertheless, recently BABAR P, 0| first observed sizable transverse 
fractions in the B (f)K* decays, where the transverse polarization amplitudes are comparable to 
the longitudinal one. This result was confirmed later by BELLE P, In other words, in terms 

of helicity amplitudes the data show that ± H p f« l-f^ooP (or 2|A||p 2|^_lP ?s |^oP 

in the transversity basis). Such an anomaly in transverse fractions is rather unexpected within 
the SM framework. Efforts have already been made for finding a possible explanation in the SM 
or NP scenario. In the SM, according to Kagan [3], the nonfactorizable contributions due to 
the annihilation could give rise to the following logarithmic divergent contributions to the helic- 
ity amplitudes: Hoo,H— ~ O [{I / ml) In^ {mi, /At)], H++ ~ 0[{l/m^)\n^{mb/Ah)], where A,, is 
the typical hadronic scale. This in turn may enhance the transverse amplitudes required to ex- 
plain the anomaly. However, in the perturbative QCD (PQCD) framework, Li and Mishima 0] 
have shown that the annihilations are still not sufficient to enhance transverse fractions. An- 
other possibility for explaining the polarization anomaly advocated by Colangelo et al. Q is the 
existence of large charming penguin and final state interaction (FSI) effects. However they got 
|AoP(i? — > pK*) < \Ao\'^{B 4>K*), in contrast to the observations lal, where the 

normalization J2i l^iP = 1 is adopted. With the similar FSI scenario, Cheng et al. |lO| obtained 
l^oP : l^iiP : = 0.43 : 0.54 : 0.03, which is also in contrast to the recent data^3,0|- Now 

the question is: Is it possible to explain this anomaly by the NP? If yes, what types of NP oper- 
ators one should consider? Some NP related models have been proposed where the so-called 
right-handed currents 57^ (1 + 75)6 87^^(1 ±75)3 were emphasized j7|. If the right-handed currents 
contribute constructively to A± but destructively to Aq m, then one may have larger |A_|_/Ao|^ to 
account for the data. However the resulting \A^^\- « |lp i will be in contrast to the recent 
observations 3,0]. See also the detailed discussions in Sec. Hll 

In the present study, we consider general cases of 4-quark operators. Taking into account all 
possible color and Lorentz structures, totally there are 20 NP four-quark operators which do not 
appear in the SM effective Hamiltonian (see Eqs. H30() and H31|) ). After analyzing the helicity 
properties of quarks arising from various four-quark operators, we find that only two classes of 



2 



four-quark operators are relevant in resolving the transverse anomaly. The first class is made of 
operators with structures cr(l — 75) — 75) and (1 — 75) (8) (1 — 75), which contribute to different 

helicity amplitudes as Hqq : H : ~ ©(l/mf,) : 0{l/m1) : 0(1). The second class consists 

of operators with structures a{l + 75) ® (7(1 + 75) and (1 + 75) ^ (1 + 75), from which the resulting 

amplitudes read as Hqq : : H ~ 0{l/mb) : 0{l/m1) : 0(1). Moreover, the above (pseudo- 

)scalar operators can be written in terms of their companions, the (axial-)tensor operators, by Fierz 
transformation. Finally, there is only one effective coefficient relevant for each class. We find that 
these two classes can separately satisfy the two possible solutions for polarization phase data, which 
is due to the phase ambiguity in the measurement, and the anomaly for large transverse fractions 
can thus be resolved. The tensor operator effects were first noticed by Kagan 0] (see Sec. O for 
further discussions). 

The organization of the paper is as follows. In Sec. ^ we first introduce the SM results for the 
polarization amplitudes in the B — > (f)K decay within the QCD factorization (QCDF) framework. 
After that we give a detailed discussion about how the NP can play a crucial role in resolving the 
large transverse polarization anomaly as observed by BELLE and BABAR. The reason for choosing 
the two classes of operators with structures (i) a{l — 75) cj(1 — 75), (1 — 75) ® (1 — 75) and (ii) 
cr(l +75) (8>cr(l -|-75),(1 +75) (8> (1 +75) is explained and the relevant calculations arising from 
these operators are performed. We discuss the possibility for the existence of right-handed currents 
57^(1-1-75)6 S7'^(lib75)s which was emphasized in From the point of view of helicity conservation 
in the strong interactions, we discuss various contributions originating from the chromomagnetic 
dipole operator, charming penguin mechanism, and annihilations. Some observables relevant in 
our numerical analysis are defined in this section. In Sec. II I II we summarize input parameters e.g. 
Kobayashi-Maskawa (KM) elements, form factors, meson decay constants, required for our study. 
Sec. IIVI is fully devoted to the numerical analysis. We discuss in detail two scenarios, which are 
separately consistent with the two possible polarization phase solutions in data due to the phase 
ambiguity. We obtain the best fit values for the NP parameters which can resolve the polarization 
anomaly. Numerical results for observables are collected in this section. Finally, in Sec. we 
summarize our results and make our conclusion. 



II. FRAMEWORK 



A. The Standard Model results in the QCD factorization approach 



The best starting point for describing nonleptonic charmless B decays is to write down first the 
effective Hamiltonian describing the processes. The processes of our concern are the B — > (j)K* 
decays which are penguin dominated. In the SM, the relevant effective weak Hamiltonian Ti^.^ for 
the above AS = 1 transitions is 



7^ -^^ 



(1) 



Here Cj's are the Wilson coefficients and the 4-quarks current-current, penguin and chromomagnetic 
dipole operators are defined by 
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current-current operators: 

O" = {ub)v-A{su)v^A 
01 = {ch)v-A{'sc)v-A 

QCD-penguin operators: 

O3 = {sb)v-AY^{qq)v-A, 

O5 = ish)v^AY^{qq)v+A, 
1 

electroweak-penguin operators: 

3 

O7 = ■::^b)v-A^eq{qq)v+A, 



O2 = {Uabi3)v-A{sf3Ua)v-A, 
O2 = (c«6/3)y-A(S/3Ca)\/-A, 



O4 



{sabp)v-AY^{qpqa, 



V-A, 



{Sab(})v^A^{qpqa)v+A, 



Os = ■::{sabp)v-A^eq{qpq, 



a)V+A, 



O9 = ^{sb)v-A^eg{qq)v-A, Oio = -{Sabf3)v-A^eg{q0, 



(2) 



(3) 



a)V-A, 



(4) 



chromomagnetic dipole operator: 



(5) 



where a, (3 are the SU{3) color indices, V ± A correspond to 7^(1 ± 7^), the Wilson coefficients 
Cj's are evaluated at the scale //, e and g are respectively QED and QCD coupling constants and 
T"'s are SU{3) color matrices. For the penguin operators, O3, . . . , Oio, the sum over q runs over 
different quark flavors, active at /x ~ mb, i.e. q €{u, d, s, c, b}. 

In the present work, we will embark on the study of -B — > (pK* decays in the approach of the 

ly a] 

Gf 



Q I 

QCDF. The B (pK decay amplitude with the (p meson being factorized |1J| reads 



V2 



i-vtbVt:) 



as + 04 + ^5 - ^("7 + 09 + aio] 



X 



(6) 



which is penguin dominated. The annihilation contribution which is power suppressed is neglected 



here The (pK*^ decay amplitude can be obtained by considering CP transformation. 

As far as the charged i?-meson decay is concerned, the dominant contribution also comes from 
the penguin operators, while the contribution due to Oi,02 is color and KM suppressed. In the 
scenario, where (j) is factorized, the decay amplitudes for B^, B^, B^, B^ are almost the same. The 
factor ^' ''^) in Eq. © is equal to 

x(b"k'','I>) = {^(q,e,)\{ss)v-A\0){K*\p',e2)\isb)v-A\B\p)), 

"2^ *^J.,*|y„a/^3^r^„2', 



ms + rriK 



-e,uapere*2''pyV{q' 



{mB + mK*)el ■ £2^1(9^) - (e^ •p)(e2 ■ p) 



2A2{q^ 



niB + rriK-* 



(7) 
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where the decay constants and form factors are defined by 



p.* 



{K*\p\e2)mB\p)) 



{K*\p\e2)W\B\p)) 



{tub + mK*)e2^Ai{q ) - {el ■ p){p + p) 



-2imK- 



(8) 



with niB and rriK* being the masses of B and K mesons, respectively, q = p — p', ^43(0) = Ao(0), 
and 



Asiq') 



ruB + rriK* ^ 2^ nT-B - niK* . , 2n 



-AM) 



2mK* 2mK* 
It is straightforward to write down the decay width, 



-A2iq' 



Pc 



Sirm 



B 



ool^ + \B++\'^ + \H- 



-*o 



(9) 



(10) 



where Pc is the center mass momentum of the (j) ox K meson in the B rest frame. Hqq^ 
H are the decay amphtudes in the hehcity basis and in QCDF ^, they are given by 



Gp 



Hoo = —^{Vtbyt*s)asM{^U'm'<t>){mB + mK*) \aAi{m^) - bA2{m^ 



Gf. 



^ niB + nriK* 



(11) 



with the constants a = (m^ — ""t-^ — '"T'i-*)/(2m0mx*), b = {2'm^pf.)/[m^mx*{mB ]■ Here 

'^SM = 03 + 04 + fls — ^{0.7 + Og + Oio)- The superscript h in a^'s denotes the polarization of (j) and 
— *o 

K mesons; /i = is for the helicity 00 state and h = ± for helicity ztit states. Note that the weak 
phase effect is tiny in a^j^ and is thus neglected in the study. Such helicity dependent effective 
coefficients cigj^^j do arise in the QCDF, however in the naive factorization (NF), they turns out to 
be same, i.e. a^j^ = a^j^ = a^j^ = asu- In the NF, one can rewrite the above amplitudes in the 
transversity basis as 



A 



■SM 



Gp 
72 



{VtbVts)asM{if4,'m^){mB + rriK*) aAi(m^) - bA2{m^^] 



We choose the coordinate systems in the Jackson convention, consistent with what BaBar and Belle 
did Q|. In the B rest frame, if the z axis of the coordinate system is along the the direction of the flight 
of the 4> meson and the transverse polarization vectors of </> are chosen to be e^(±l) = (0, =Fl, — i, 0)/\/2, 
then the transverse polarization vectors of K* are given by e^. (±1) = (0, ^pl, +*, 0)/\/2 in the Jackson 
convention, but become e^. (±1) ~ (0, ±1, — i, 0) / \/2 in the Jacob- Wick convention. Therefore in the NF, 
arg(A|| _l/Ao) equal to tt in the Jackson convention, but are zero in the Jacob- Wick convention. Note that 
in the two conventions, the longitudinal polarization vectors are the same as €^(0) = (pc, 0, 0, E^) / and 
e^. (0) — {pc, 0,0, —EK')/mK*- Here the amplitudes satisfy A = Aq -|- Ay +A±,A = Aq+A^^ — A±, where 
the kinematic factors are not shown. 
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^il'^ = -^{VtbVt*)asM{i\^Um^){mB + 'mK*)Ai{rnl), 
In the QCDF, a^'s are given by 

+(C4 + ce) ^ (G'^(.i) + lies + cio) E (G'^(^i) - ^ 



07 



^ -L '^8 Ois Cf (i^h , fh ^ Mnh 



h , C7 

as = ^8 + 7^' 

where Cp = {N^ — l)/{2Nc), Si = mf/ml, Xq = VqbV*^,, and q' = d, s. Note that we have given the 
expressions for oj", which may be relevant for charged B decays, arising due to Oi and O2 in 
Eq. (^. There are QCD and electroweak penguin- type diagrams induced by the 4-quark operators 
Oi for i = 1,3,4,6,8,9,10. These corrections are described by the penguin-loop function G^{s) 
given by 

G°(s) = In— + 4/ du^fju) dxx{l-x)\n[s-ux{l-x)], 

3 3 mi) Jo " Jo 

G±(.) = ?-?lnA 

3 3 m-b 

+2 1 du[ gl^""' (u) ± 7^^V^^ dx x{l - x) ln[s - ux{l - x)] . (14) 



du 

In Eq. 1)13(1 we have also included the leading electroweak penguin-type diagrams induced by the 
operators Oi and O2, 



C = ( fG^isu) + fG^'isc)) + (15) 
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The dipole operator Os„ will give a tree-level contribution proportional to 



1 ^tiu 
du 







l-u ' 



In Eq. (|13|) . the vertex correction is given by 



du 4 



= -121n-^ - 18 + /^ 

where we have used the naive dimensional regularization (NDR) scheme [if 

li.lulxi'^ - 75) ® 7^7''7^(1 - 75) = 4(4 - e)7^(l - 75) ® 7^(1 - 75), 
7m7/^7a(1 - 75) ® 7^(1 - 75)7''7'' = 4(1 - 2e)7/,(l - 75) ® - 75), 
7/.7^7a(1 - 75) ® 7^7V(1 + 75) = 4(1 + 5)7^(1 " 75) ® 7^(1 + 75), 
7/.7^7a(1 - 75) ® 7^(1 + 75)7''7^ = 4(4 - 4e)7^(l - 75) ® ^^"{1 + 75), 



(16) 



(17) 



(18) 
(19) 
(20) 
(21) 



with D = 4: — 2e, and have adopted the MS substraction. An explicit calculation for arising 
from vertex corrections, yields 



ft 







1 , / 1 — 2x 
dx ^fAx) 3— Inx — Sin ] , 



I — X 



dx [ g'i^''\x) ± 



4 dx 



3 In x — Sivr ) . 

1 — X 



(22) 



The hard kernel f^j for hard spectator interactions, arising from the hard spectator interactions 
with a hard gluon exchange between the emitted vector meson and the spectator quark of the B 
meson, have the expressions: 



III 







dp 



di] ■ 



P Jo 



dr] 



di 



47r2 2fBfK*rnK* 



rin?^h± 







p Jo 



K*ia) 



dr] 



di 



where 



+ ^J^^ dv{^\{v) - gf\v))V 



(23) 



2mBPc yBK 



rriB + rriK* 



mB + rriK' 
(ml). 



(24) 
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Note that can be obtained from in Eq. (|17() with the replacement of g'^^"''^ —g^^°^ ■ We 

will introduce a cutoff of order AqcD/'^^b to regulate the infrared divergence in ///. Note also 

that we have corrected ^ the QCDF results in Ref. which were done by Cheng and one of us 

(K.C.Y.). The key point for the calculation is that one needs to consider correctly the projection 

operator in the momentum space, as discussed in Appendix which may explain the difference 

with Ref. ^- In the calculation, we take the asymptotic light-cone distribution amplitudes 

(LCDAs) for the light vector mesons, and a Gaussian form for the B meson wave function. |14|. 

We now make a SM estimate for various helicity amplitudes from a power counting point of view. 

For 4>K*^ , the helicity amplitude Hqq arising from the iV — A) ®{y — A) operators is 0(1), 

since each of (j) and K mesons, with the quark and antiquark being left and right-handed helicities, 

respectively, requires no helicity flip. For H , the helicity flip for the s quark in the (j) meson is 

required, resulting in m^/rrib suppression for the amplitude. Finally in two helicity flips for s 

— — *o 

quarks are required, one in cp meson and the other in the B ^ K form factor transition, which 
cause a suppression by {m^/mi,){A/mi,), where A = uib — rrii,. In a nutshell, the three helicity 

amplitudes in the SM can be approximated as Hqq : H : -ff++ ~ 0(1) '■ 0{\/mi,) : 0{l/ml). 

One should note that for the CP-conjugated B^ — > <j)K*^ , the above result is modified to be 
ifoo : : ~ 0(1) : 0(l/mb) : 0(l/m^). The results extended to various possible NP 

operators together with the SM operators will be shown later in Table ^ a-nd also illustrated in 

Fig.m 



B. New physics: hints from the BABAR and BELLE observations 

The large transverse B (f>K* fractions as have been observed by BELLE and BABAR 0, 0, 
, may hint a departure from the SM expectation for the longitudinal one. Within the SM, the 



QCDF calculation [ij] yields 

l-Ro = 0{l/ml) = Rt , (25) 

where = \'Ao\'^ /\Atot\'^ and Rt = R\\ + R± = (P||P + |A^P)/ptotP- The observation of large 
Rt, as large as 50%, may be possible to be accounted for in the SM, but here we are considering 
the new physics alternatively ^. In the transversity basis the recent experiments 0, 0] have shown 
that 

\Ao\\= \Hoof) ^ \ATf{= Pill' + \A^f) , (26) 



^ We especially thank A. Kagan for pointing out that some terms in Cg may be missed in Ref. We 

were therefore motivated to recalculate the QCDF decay amplitudes. 
^ Our G+ does not agree with that obtained by Yang et al. . ^ 
^ In the PQCD approach, annihilation contributions appear to be too small to resolve the puzzle 8], but 
Li has recently argued that an decrease in one of the B K* form factors could be helpful. Nev- 
ertheless, using the QCDF, Kagan showed that the suppressed annihilations could account for the 
observations with modest values for the BBNS parameter pA- See also the discussion after Eq. H35() . 
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wherein = + i/__)/\/2 and = - i/— )/\/2. One may need NP to explain such 

a large Rt (~ Rq)- A set of NP operators contributing to the different helicity amplitudes like 

Hoo : H.. : H++ ~ Oil/mt) : : 0(l/mg) , (27) 

or 

Hoo ■■ H— : H++ ~ : 0(l/mf ) : , (28) 

could resolve such polarization anomaly. Note that H (in the former case) and (in the 

latter case) are of 0(1), while Hoo is always of 0{l/miy), in contrast to the SM expectation. The 
detailed reason will be seen below. 



1. New physics operators 

Now Eqs. (|27j) and (|28j) can serve as guideline in selecting NP operators. To begin with, consider 
the following effective Hamiltonian TL^^: 

p 30 

n''^ = ^Y.i^ii^')o,{^^)] + H.c., (29) 

which may be generated from some NP sources and contains the following general NP four- 
operators: 

• four-quark operators with vector and axial-vector structures: 

On = 57^^(1 + 7^)6 S7^(l + 7^)5 , Ou = Sa^f^il + i^)bp 5^37^(1 + -i^)so, , 
Oi3 = 57^^(1 + 7')^ S7/.(1 - l')s , Ou = Sc.j'^il + l'')hp sp^^H - 7')s„ , 

(30) 

• four-quark operators with scalar and pseudo-scalar structures: 

Ol5 = s{l + j^)b s(l + 7^)5 , O16 = s«(l + 7^)&/3 s/sil + 7^)Sa , 

On = s{l - 7^)6 s(l - 7^)5 , Oig = s«(l - 7^)6/35 sp{l - j^)so, , 

Oi9 = s(l + 7^)b s(l - , O20 = s«(l + 7^)^/3 s/3(l - 7^)Sa , 

O2I = S{1 - 7')ft S(l + 7')^ , O22 = S«(l - 7')^'/3 S/3(l + 7')Sa , (31) 

• four-quark operators with tensor and axial-tensor structures: 

O23 = saf^'il + 7')^ sa^,{l + 7^)5 , O24 = Saa^'il + 7')^'/3 S/3a^.(l + 7')sa , 
O25 = sa'^"(l - l^)b sa^uil - 7')^ , O26 = s<,a'^'^(l - 7')fc/3 S/3CT^.(1 - 7')5a , 
O27 = sct'^"(1 + 7')fe sa^,{l - 7^)5 , O28 = s„a^'^(l + 7^)6^ s^a^,(l - 7^)5^ , 
O29 = sa^^il - l^)h sa^^(l + l^)s , O30 = s„a^'^(l - 7^)6^ S/3CT^^(1 + 7^)sa • 

(32) 
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Here Cj with i = 11, ... ,30 are the Wilson coefficients of the corresponding NP operators and /.i 
the renormalization scale, chosen to be mi, here. Now we give an estimation of several types of NP 
operators, contributing to various B — > (pK helicity amplitudes. In Fig. [TJ we draw the diagrams 
in the B rest frame, where qi, are the s quarks, and q2 is the s quark, (^i, ^2) and (^3, ^4) form cp 
and K* , respectively, q^ is the spectator light quark which has not any preferable direction. qi,q2, qs 
are originated from the following operators: gsFii qiT2q2 for Os-g, O11-14, O23-30 or qiTib q3T2q2 
for Oi5-22- If the helicity for qi or q2 is flipped, then the amplitude is suppressed by a factor of 
iT^cfi/iT^b- On the other hand, if the helicity of q^ is further flipped, the amplitude will be suppressed 
by {mfj,/mb){A/mi,), with A = — m^. The results are summarized in Table 



q4 b _« qj 

" — o :_ 



q4 b ^ qi 

^ n 

'i3'^ q2 



q4 



b 



q3 



o, 



q4 



b 



13 



o, 



q4 



^3 



b 



o 



15,16,23,24 



q4^ 

qj 



b 



^li 

^2 



q4 



^3 



b 



o 



19,20,27,28 



^^1 



q4^ 

qa 



b 



1- ,qi 

O "^2 

^21,22,29,30 



FIG. 1: The main helicity directions of quarks and antiquarks arising from various four-quark 
operators during the B decay, where the solid circle denotes the b quark, qi,q3 are the s quarks, 
and q2 is the s quark. ((71,92) and ((73,^4) form 4> and K* , respectively. 94 is the spectator light 
quark which has no preferable direction. The short arrows denote the helicities of quarks and 
antiquarks. See the text for the detailed discussions. 

From the Tabled we see that both (pseudo-) scalar operators Ois-is and (axial-)tensor operators 
O23-26 satisfy the anomaly resolution criteria as given by Eqs. (|27j) and (|28jl . while the rest are 
not. However, through the Fierz transformation, it can be shown that Ois^ie and Oi7,i8 operators 
can be expressed as a linear combination of 023,24 and 025,26 operators, respectively, i.e., 

O15 = Y^023 — g024, 
^16 = 7^024 — c^23! 
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TABLE I: Possible NP operators and their candidacy in satisfying the anomaly resolution criteria. 
We have adopted the convention Pi (g) r2 = sPi^ ST2S. 



Model 


Operators 








Choice 


SM 


7/^(1 -75) 07^(1^75) 


0{l) 


o(iM) 


0{l/ml) 




NP 


7''(l + 75)«'7/.(l + 75) 


0{l) 


0{l/ml) 


0{l/mb) 


N 


NP 


7''(1 + 75) ^7^(1-75) 


0{l) 


0{l/ml) 


0(1M) 


N 


NP 


(1 + 75)0(1+75) 


0{l/m,) 


0{l) 


0{l/ml) 


Y 


MP 


iJ- — 75J U — 75j 






n(^ \ 


V 

1 


NP 


(1 + 75) 0(1- 75) 


0{1) 


0{l/ml) 


0(1M) 


N 


NP 


(1-75)0(1+75) 


0{1) 


0{l/m^) 


0{l/ml) 


N 


NP 


^^''(1+ 75)0^^^(1 + 75) 


0{l/m,) 


0{l) 


Oil/ml) 


Y 


NP 


^^"(l-75)®<T^.(l-75) 




0{l/ml) 


0{l) 


Y 


NP 


fT^''(l+75)«'<^Mi^(l-75) 


0{1) 


0{l/ml) 


0(1M) 


N 


NP 


(7^''(1- 75)0^^^(1 + 75) 


0{1) 


0(1M) 


0{l/ml) 


N 



On — Y^C)25 - g^26, 

O18 = ^02G-\025. (33) 

Before we continue the study, five remarks are in order, (i) The sa^'^{l + 75)6 scj^,y(l + 75)5 
operator, which could maintain « ^"cst mentioned by Kagan [7| ^. (ii) We do not 

consider NP of left-handed currents, 57^(1 — 75)6 S7'^(l + 75)s, which give corrections to SM Wilson 
coefficients, ci_io, since they have no help for understanding large polarized amplitudes and are 
strongly constrained by other B PP,VP observations (iii) On -14 are the so-called right- 
handed currents, emphasized recently by Kagan These operators give corrections to amplitudes 
as 

^0,11 {-VtbVtDasM - {ail + ai2 + aia), 

A± oc {-VtbVt*s)asM + (an + ai2 + ais), (34) 

where 

ail = cii + ^ + nonfact., 012 = C12 + ^ + nonfact., 013 = C13 + ^ + nonfact., (35) 

with "nonfact. "= nonfactorizable corrections. Note that 011,12,13 enter the ^o,|| amplitudes with 
a "minus" sign due to the relative sign changed for Ai,A2 form factors as compared to the SM 
amplitudes in Eq. (|llj) . If the right-handed currents contribute constructively to A± but destruc- 
tively to ^0,11 > then one may have larger I^^/AqP to account for the data. According to the 
SM result |]4_l/AoP ~ 0.02 in Eq. (jSSI), we need to have |aii + 012 + ai3|/| - VtbVtlasM\ ~ 1.5 



^ However, the contributions arising from the sa^'^{l + 75)6 sa^^{l —75)5 operator to different polarization 
amplitudes should be Hqq : H-- : H++ ~ 0{1) : 0{l/ml) : 0{l/mb), not as mentioned in |2l- 
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such that |^j_/^oP — 0.25. However the resulting l^yP <C will be in contrast to the 

recent observations 0, 0]. (iv) Since, in large limit, the strong interaction conserves the 
helicity of a produced light quark pair, helicity conservation requires that the outgoing s and 
s arising from s — s — n gluons vertex have opposite helicities. The contribution of the chro- 
momagnetic dipole operator to the transversely polarized amplitudes should be suppressed as 

Hqo : H : ~ 0(1) : 0{l/mh) : 0{l/ml); otherwise the results will violate the angular 

momentum conservation. Actually if only considering the two parton scenario for the meson, the 
contributions of the chromomagnetic dipole operator to the transversely polarized amplitudes equal 
to zero ^ (see Eq. (jl6p and Appendix^ for the detailed discussions). Similarly, the s and s quark 
pairs generated from c, c annihilation in the charming penguin always have opposite helicities due 
also to the helicity conservation. Hence, that contributions to the transversely polarized amplitudes 
are relative suppressed, in contrast with the results in Refs. (v) With the same reason as the 

above discussion, in the SM, the transversely polarized amplitudes originating from annihilations 
are subjected to helicity suppression. A suggestion pointed out by Kagan 0] for the polarization 
anomaly is the annihilation via the {S — P) <Si {S + P) operator, which contributes to helicity 

amplitudes as Hqo,!! ~ 0[{l/ml)\n^{mb/Ah)], 11++ ~ 0[{l/ml)\v?{mi,/ Kh)]. However this 

contribution to H is already of order l/mf,^ although it is logarithmic divergent. 

We now calculate the decay amplitudes for B —>■ 4>K due to Oi5_i8 and O23-26 operators in 
Eqs. (EH) and (EH). The amplitudes for 5° cpK*^ can be obtained by CP-transformation. The 
matrix elements for (axial-) tensor operators 023,25 can be recast into 

{ct>{q,ei)X{p\e2)\-sa^''{l±l^)s sa^^il ±-f^)h\B{p)) =(l+ ^ 



2N, 

xfJ^Se^^paeTe^pPp"^ T^{q^) ^ AiUq^) [{el ■ e^)(m| - m|,.) " 2(6^ • p) (e^ • p) 

±8iUq^){el-pM-p) .""K ] , (36) 

under factorization, where the tensor decay constant /J is defined by 2^, 23, 24| 

(0((?,6i)|sa^'^s|O) = -iflieTq" - e^q^), (37) 

and 

{K\p',e2)\sa^,q-{l+^^)b\B{p)) 

= ie,up.e*W 2ri(g2) + T2{q') {el^{ml - mj,,) - {e* ■ p) {p + p'),} 

+ Uq'M-pB) jg^ - 2 2 iP + P)A , (38) 

with 

ri(o) = r2(o). (39) 

The helicity amplitudes for the B^ decay due to the NP operators are (in units of Gf/\^) given 
by 

^00^ = --^ifl^B (^23 - 025) \h2T2{ml) - hsT^iml)] , 
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H±± = -^ifjml {d23 [±/iri(m^) - /2r2(m^)] + 025 [±/iTi(m^) + /2r2(m^)] } 



(40) 



and in the transversity basis, the amphtudes becomes (in units of Gp/V^) 

= -AifJniB [023 - 025] /i2?2(m^) - hsT3{ml) 
a'''' = 4i^/2/Jm|(a23-a25)/2T2(mi), 



A 



NP 



4^\/2/Jm|(a23 + a25)/iTi(m? 



(41) 



where 



2pc 
rriB 



, /2 



m 



B 



m 



K* 



m 



B 



1 



and 



^23 



025 



1 + 
1 + 



1 



2N, 
1 



1 

1 

C23 + Y^cis 




m 



B 




m 



(42) 



K* 



1 



6 

1 1 

C25 + TTTCit - -C18 



C16 + 



12 



6 



+ 



1 1 
1 1 

iv: + 2 



C24 + Y^cie 



C15 I + nonfact., 



C26 + ^ci8 - ^ci7 ) + nonfact., 



(43) 



are NP effective coefficients defined by 023 



1 023 1 e*''23e»</'23 



',025 = I025I 6**^25 e*'P25 with 4>23,4>2!i being 
the corresponding NP weak phases, while 523)<^25 the strong phases. Note that here we do not 
distinguish effective coefficients for different hehcity amphtudes since those differences are relatively 
tiny compared with the hierarchy results in Eqs. (|27() and H28|) . A further model calculation 
for 023,025 will be published elsewhere |25|. Note that if applying the equation of motion to 4- 
quark operators in deriving the matrix in Eq. H41|) . we can obtain the following relations: Ti ~ 
VmB/imB+mK*),T2 c:i_A- ^ m q /_(rnB — mK*),T?, ~ A2, consistent with results by the light-cone sum 
rule (LCSR) calculation 
the results of the — > 



Ai ^ m n l (m 
3,013. 



The — > (j)K*^ polarization amplitudes can be obtained from 
decay by performing the relevant changes under CP-transformation. 
The total SM and NP contributions for the B^ and decays can be written as 



A{BO 
A{B^ 



*0x 



(j)K 



A{BO 
A{B^ 



(pK 



*0^ 



)SM 



+ A{B^ 



*0\ 



)SM 



4>K*^)np, 
+ A{B^ ^<t>K*'^)NP. 



(44) 



With these decay amplitudes in the transversity basis, we can evaluate physical observables: 



l^olM^iilM^j 



l^olM^ii 



^00, A_L_|_, A^o, ^±11, ^110, ^00, i;_L_L, i;_L0, Sii 



and the triple products ^^n, AJj,, Aj^,Arp |2a]. The observables Ahh and Y^hh are defined as 

Ahh = \ (lAP + \Ahf) , 
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A±i = -Im{A±A; -A^_A*), 
A||o = ReiAiiAl + A^Al), 

= -Im{A±A* +A^_A*i), 
S||o = Re{AiiA*o-A\\Al), 

with h = 0, ||,_L and i = 0, ||. Here we adopt the normalization conditions X^J^i 
X]j l^jp = 1- The two triple products Aj^ and A^)p are defined as 

^0 ^ ImiA^A*^) 



(45) 
1 and 



An 



A< 



Im{A_i_Al) 



\Ao\'+ A 



(46) 



In our numerical analysis, we will focus on the studies of these quantities. The CP-conjugated 
and A^' can be obtained by replacing ^0)^||)^± by their CP-tranformed forms ^0)^||)^±- 
Observables like '^w, Syo, A_|_j (with A = 0, ||, _L and z = 0, ||) are sensitive to the NP 2^], which, in 
absence of the NP, strictly equal to zero. The triple product Aj, or A^^ can exhibit the relative phase 
between A± and Aq or between A± and j4||. The differences between Aj, and their CP-conjugated 
parts, i.e. A_|_j = — (with i = 0, ||), are CP-violating (and also T-violating following from the 
CPT invariance theorem) quantities. Therefore, any non-zero prediction of Syg, A^j resembles 
the evidence of a new source of CP-violation. Moreover, since CP-violated effects are expected to 
be negligible within the SM, sizable A_|_o or A^y may also imply the existence of the NP. We will 
look for these possibilities from a detailed numerical study. 

III. INPUT PARAMETERS 

The decay amplitudes depend on the effective coefficients Oj's, KM matrix elements, several 
form factors, decay constants. 

A. KM matrix elements 



We will adopt the Wolfenstein parametrization, with parameters A, A, p and rj, of the KM matrix 
as below 

f Vud Vus Kb \ / 1 - \>? A A\\p - ir,) \ 

Vcd Vcs Vcb = -A l-iA^ AA2 



Vkm 



\Vtd Vt 



ts 



Vt 



tb J 



A\^{1 - p-iri) -A\^ 



J 



We employ A and A = sin^c at the values of 0.815 and 0.2205, respectively, in our analysis. The 
other parameters are found to be p = p(l - = 0.20 ± 0.09 and 77 = 77(1 - /2) = 0.33 ± 0.05 
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B. Effective coefficients Oj, form factors and decay constants 



The numerical values for the effective coefficients with i = 1,2, 10, which are obtained in 
the QCDF analysis Q], are cataloged in Table ITII The effective coefficients are the same both 



for B and B, but not so for a. 



(-,+) 



In the third and fifth columns of Table ITTl the a^'s with the 



superscript being bracketed are for the B —>■ (j)K process, otherwise for the B (j)K* process. 



TABLE II: Effective coefficients for B 
analysis 3|) where "'^''^ ^ 



''2,3 



(f>K*{B — > (f>K*) obtained in the QCD factorization 
sensitive to the nonfactorizable contribution /^^ which is opposite in 



sign to fjj , are obviously different from 02 3, a2 



,(+) 



a? 


1.0370 + 0.0135i 




1.0900 + 0.0187f 




1.0180 + 0.0135i 


4 


0.0764 - 0.0793i 




-0.2351 - 0.1096i 


02 


0.1898 - 0.0793i 




0.0055 + 0.0026i 




0.0156 + 0.0035i 




0.0019 + 0.0026i 


al 


-0.0347 - 0.0068i 




-0.0366 + O.OOli 


04 


-0.0310 - 0.0036i 


«g 


-0.0046 - 0.0030i 




0.0023 - 0.0019i 


-(+) 


-0.0077 - 0.0030i 


a? 


-0.0001 - O.OOOli 




0.00001 - O.OOOH 


-(+) 


0.00009 - 0.00005f 




-0.0092 - 0.0003i 




-0.0096 - 0.0002i 


-(+) 


-0.0088 - 0.0002i 


"10 


-0.0004 + 0.0006i 




0.0023 + O.OOlOi 




-0.0014 + 0.0007i 



For the decay constants, we use |1 



J 6 



237 MeV, fx* 



160 MeV, and 



Jb = 190 MeV. For the B K* transition form factors, we adopt the LCSR results in j2^ with 
the parametrization 



F{q')=F{0) exp ci(gVm^) + C2((?V"^i 



2 ^2 



(47) 



which were rescaled to account for the B K*j data. The values of the relevant form factors 
and parameters are given in Table ITTTI The reason for choosing this set of form factors is because 



TABLE III: The values for the parametrization of the B 
renormalization scale for Ti, T2, T3 is = m^. 



K* form factors in Eq. ^ The 





yli(O) 


^2(0) 


^o(O) 


y(o) 


ri(o) 


^2(0) 


^3(0) 


m 


0.294 


0.246 


0.412 


0.399 


0.334 


0.334 


0.234 


Cl 


0.656 


1.237 


1.543 


1.537 


1.575 


0.562 


1.230 


C2 


0.456 


0.822 


0.954 


1.123 


1.140 


0.481 


1.089 



the Ti (0) value extracted from the B K*^ and B Xg'j data seems to prefer a smaller 
one 
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IV. NUMERICAL ANALYSIS 



We will estimate the NP parameters which may resolve the polarization anomaly in B{B) 



(j)K*[K ) decays An enhancement in transversely polarized amplitudes by 50% can therefore 

take place in our NP scenario since the SM polarization amplitudes, Aq : Ay : Aj_ ~ 0(1) : 

/VP IVP IVP 

0{l/mb) : 0{l/mb), are modified to be Aq : A\^ : A^ ~ 0{l/mb) : 0(1) : 0(1), as given in 
Eq. (|lT|l which allows us to find solutions in the NP parameter space (|a23| , 623, 4>23, |a25| , ^25, 4'25) 
for explaining the B — > (f)K* polarization anomaly. 

Choosing the normalization conditions J2i l^jP = J2i l^iP = Ij ^^'^ setting arg(^o) = 
arg(^o) = 0, one can measure the magnitudes and relative phases of the six ^o,||,±) ^o,||,± PO" 
larization amplitudes ^, giving 8 measurements, and then extracts 12 observables in Eq. H45|) as 
well as the triple products in Eq. (|46|) . We take the average of the B AB AR and BELLE data in our 

analyses for estimating NP parameters and consequently obtain the predictions for observables. 
For simplicity, we neglect the correlations among the data. The xf of any observable Oj with the 
measurement Oi{expt) it {laijexpt is defined as 

2 

(48) 



Oi{expt) — Oiitheory) 



For different observables, the total equals to x^ = Y^f=iXi.- ^ \i&'sX, fit analysis, we 

consider the following 8 observables: 

|Aop,|^±P,|:4op,|:4x|2,arg(^l|),arg(Ax),arg(:4||),arg(Ax), (49) 

as our inputs. For the purpose of performing the numerical analysis easily, we have converted the 
BABAR measurements into the above quantities, as shown in Tables [TVl and IVl 

Since the interference terms in the angular distribution analysis j^, |^ are limited to i?e(A||Ao), 
/m(j4^^g), and Imi^A^jC^^ there exists a phase ambiguity: 

arg(^ll) ^ -arg(]4||), 

arg(^_L) ±vr - arg(^_L), 
arg(:4x) - arg(:4||) ^ ivr - (arg(:4^) - arg(:4||)). (50) 

Therefore, the world averages for arg(A||) and arg(y4_L), given in Tables HVl and Ivl can be 

arg(:4||) = -2.33 ± 0.22, arg(A_L) = 0.59 ± 0.24, (51) 

or, following from Eq. (|50() . 

arg(A||) = 2.33 ± 0.22, arg(A_L) = 2.55 ± 0.24. (52) 

From Eq. H51() . the phase difference for Ax_ and A\^ reads 

arg(A±) - arg(^ii) « tt, (53) 



^ For simplicity, in the present study we have chosen the convention arg(ylo) = arg(Ao) = 0, i.e. we do not 
consider here the physics arising from the difference between arg(^o) and arg(Ao). 
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but, on the other hand, from Eq. ()52() . becomes 

arg(^±)-arg(A||) RiO. (54) 
The resultant imphcations in Eqs. H53|) and 1)54(1 are discussed below. Numerically, SM and NP 

Q 

amplitudes in the transversity basis for B (pK are given by 

= -0.00307 - 0.00074 i, 
A^^^ = +0.00048 - 0.000064 i, 

A^^ = +0.00040 - 0.000073 i, (55) 



and 



aH^ ~ 2.2(|523|e*('^23+</'23) _ |^25|e^('525+te)) 





A^^ ~ -11.9 (|523|e^('^23+</'23) _ |a25|e^('^25+</'25)^ 
A^^ ~ -10.9 (|523|e^('^23+</'23) + |a25|e^(''25+</'25)^ ^ (55) 



respectively, in units of —iGp/V^- From Eq. we find that | Aj^*''^|/|^g^'^| and |^l^^|/|^o^^| are 

0.17 and 0.14, respectively. In other words, are 0{l/mh) suppressed, compared to Aq^ . On 

the other hand, the measurements for ~ l^±P ~ ^l^oPj as cataloged in Table HVl mean that 
jVP 7\rp 

A|| and A± are dominated by and Aj_ , respectively. We therefore find that the data for the 
amplitude phases in Eq. (|^T|) prefer the 025 terms in and Aj_ given in Eq. (j56|) . since there is 

a phase difference of vr between two 025 terms. Consequently, from Tabled we get 3> H if 

Oi7,i8, O25.26 NP operators are dominant. On the other hand, for the data of the amplitude phases 
in Eq. ()52() . we find that the 023 terms in Ay and in Eq. (|56() are instead favored, since 

they have the same sign. Accordingly, also from Tabled as only Oi5,i6) 023,24 NP operators are 

considered we obtain H ^ ^++1 which is consistent with the SM expectation [23,123]. Therefore 

because of the phase ambiguity, the data prefer two different types of NP scenarios: (i) the first 
scenario, where the NP is characterized by Oi7,i8,25,26 operators, while the operators Oi5^i6,23,24 
are absent, (ii) the second scenario, where the NP is dominated by Oi5,i6,23,24 operators, while 
On 18 25 26 operators are absent. 



A. The first scenario v^ith Oi5,i6,23,24 absent 



In this scenario, the NP effects characterized by Oi7,i8,25,26 operators are lumped into the single 
effective coefficient 025 = 1 025 1 e*''^^^ e^'^^s ^ where (j)2b and 825 are the NP weak and strong phases 
associated with it. Therefore, in our analysis, we have three fitted parameters, |a25|,i;^25) and 
(525. The Xmin/d-o-f- for this scenario is 4.15/5, where d.o.f. = degrees of freedom in the fit. Our 
best fit results together with the data are cataloged in Tables ITVl and IVll For illustration, we obtain 
theoretical errors by scanning the ^ Xmin + ^ parameter space. The BRs are only sensitive to 
the form factors, while the rest results depend very weakly on the theoretical input parameters and 
the cutoff that regulates the hard spectator effects in the SM calculation. To estimate the errors 
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TABLE IV: Comparison between the first NP scenario predictions and the average for BABAR 
and BELLE data yj, 0| with the phases given in Eq. dnU). The Xmin/d-o.f. for 8 inputs is 4.15/5. 



NP parameters 




NP results 


1^25 








(2.10«;1^) X 10-4 


^25 








1.15 ±0.09 


(t>25 








-0.12 ±0.09 


Observables 


BABAR 


BELLE 


Average 


NP results 


arg{A||) 


-2.61 ± 0.31 


-2.05 ± 0.31 


-2.33 ± 0.22 


-2.60 ±0.14 


arg{A||) 


-2.07 ±0.31 


-2.29 ± 0.37 


-2.16 ± 0.24 


-2.40 ±0.14 


arg(Ai) 


0.31 ±0.36 


0.81 ±0.32 


0.59 ±0.24 


0.87 ±0.12 


arg{Ax) 


1.03 ±0.36 


0.74 ±0.33 


0.87 ±0.24 


1.10 ±0.13 


|AoP,(|Ao|2) 


0.49 ±0.07 (0.55 ±0.08) 


0.59 ±0.1 (0.41 ±0.10) 


0.52 ±0.06 (0.50 ±0.06) 


0.52 ±0.04 (0.53 ±0.04) 


|A,|P,(|A|||2) 








0.20 ±0.02 (0.22 ±0.02) 


|AxP,(|Aip) 


0.20 ±0.07 (0.24 ±0.08) 


0.26 ±0.09 (0.24 ±0.10) 


0.22 ± 0.06 (0.24 ± 0.07) 


0.28 ±0.03 (0.25 ±0.02) 


a!^,(a«) 




0.28 ±0.10 (0.21 ±0.09) 




0.19 ±0.04 (0.24 ±0.03) 


Arp, (Ay) 




0.06 ±0.06 (0.04 ±0.08) 




0.09 ±0.01 (0.10 ±0.00) 



for BRs, arising from the input parameters, we ahow 10% variation in form factors and decay 
constants which may be underestimated, and the resulting errors are displayed in Table IVTl The 
NP parameters are given by 

|a25| = (2.10l[;;l^) X 10"^ = 1-15 ± 0.09, = -0.12 ± 0.09, (57) 

with the phases in radians. Note that the non-small 625 may imply that the strong phase due to 
annihilation mechanism in the SM cannot be negligible. In Tables HVl and IVTl we obtain results 
in good agreement with the data. The BABAR and BELLE [^, 0| data show that 

Aoo ^ A|||| + Ax±, 

A|||| ~ Ax±, (58) 
which can be realized as follows. The transverse amplitudes are given by 

A, = Ar + Ar, (59) 

where ^||^^ ~ —I.IA^^ in this scenario. In yly of Eq. (|59|) . the interference of A^^ and ^j^^^ is 
destructive, while in A± the interference of A^^ and A^^ becomes constructive. We thus find 
A_|_^ > A|||| and accordingly A|||| ~ 0.8A_|_^. Interestingly, because 825 + 4>25 is much closer to 
as compared to 625 — (t>2b, the above interference effects thus result in |^||p/|^_Lp < |Aj|p/|yl_Lp. 
In other words, a larger |(/)25| yields larger magnitudes of Al]jp, A^p. To get the first relation of 
Eq. ((S5|) . we first take the squares of the A\\ and A_i_ of Eq. (|59|) . and then add them up together 
with their CP-conjugated parts. The interference terms are mutually cancelled and one thus finds 
A|||| + A^_L ~ Aoo, due to \ A\\^^ ^ > l^||,± I • 

We obtain Ayg = -0.33 ± 0.04 and I1_lo = -0.44 ± 0.05, as compared with the BELLE data: 
A||o = -0.39 ± 0.14 and I;_lo = -0.49 ± 0.14. Within the SM, Ayo ~ 0.30 and S_lo ^ 0.10 are in 
contrast to the data. 
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TABLE V: Comparison between the second NP scenario predictions and the average for BABAR 
and BELLE data 0, 0| with the phases given in Eq. dnil). The xliin/d-o.f- for 8 inputs is 0.56/5. 



NP parameters 




NP results 


-|a23| 








-(1.70t«:Ji) X 10-4 


^23 = <523 — 








-0.78 ±0.10 


(t>23 








0.14 ±0.09 


Observables 


BABAR 


BELLE 


Average 


NP results 


arg{A||) 


2.61 ±0.31 


2.05 ±0.31 


2.33 ±0.22 


2.42 ±0.17 


arg{A||) 


2.07 ±0.31 


2.29 ±0.37 


2.16 ±0.24 


2.21 ±0.18 


arg(Ai) 


2.83 ±0.36 


2.33 ±0.32 


2.55 ±0.24 


2.44 ±0.17 


arg{Ax) 


2.11 ±0.36 


2.40 ±0.33 


2.27 ±0.24 


2.24 ±0.17 


|AoP,(|Ao|2) 


0.49 ±0.07 (0.55 ±0.08) 


0.59 ±0.1 (0.41 ±0.10) 


0.52 ±0.06 (0.50 ±0.06) 


0.51 ±0.04 (0.52 ±0.04) 


|A,|P,(|A|||2) 








0.26 ±0.02 (0.26 ±0.02) 


|AxP,(|Aip) 


0.20 ±0.07 (0.24 ±0.08) 


0.26 ±0.09 (0.24 ±0.10) 


0.22 ±0.06 (0.24 ±0.07) 


0.23 ±0.02 (0.23 ±0.02) 


A^,(AO) 




0.28 ±0.10 (0.21 ±0.09) 




0.22 ±0.04 (0.27 ±0.04) 


Arp, (Ay) 




0.06 ±0.06 (0.04 ±0.08) 




o.oit«:;;; (o.oiiH?) 



The consistency between data and this NP scenario requires the presence of a large strong phase 
^25 and a (smah) weak phase (f)2b- Our numerical predictions for the rest NP related observables 
are A_|_||, Sqo, Syn, ~ ±(1 — 2)%, which are marginal sensitive to (/>25. Since our analysis 
yields ^if- ~ 0.09 and yl|, ~ 0.10, we therefore obtain S^y = —A}\^ — ~ —0.19 and A_|_|| = 
Ay — Ay ~ —0.01. We observe that A_|_o^_l|| , Acp(i3 — > (/liiT*), have the same sign as 
025, whereas Soo,||||,||05 opposite sign. For a small (^25 = —0.12 (= —7°), we get 

NP related quantities A_|_o ~ ~^||o ~ ^cp ^ —0.05 which may become visible in the future 
once the experimental errors go down. It is interesting to note that the existence of the non- 
zero NP weak phase (j)25 may be hinted by the BABAR measurements of arg(A_|_ — Ax) 7^ and 
arg(A|| — A||) 7^ 0. If taking alone the BABAR data as inputs, we obtain = (16ib7)°, which could 
cause sizable effects in observations: S||o(= Re{A\\Al -~A\^l)) = 0.15 ±0.07, A9p{= 0.26 ±0.03) ^ 
A^(= 0.11 ± 0.06), S_Lo(= -A^ - A^) = -0.38 ± 0.08, A_lo(= Ay - A%) = -0.15 ± 0.07 and 
Ajlp ~ -2A^p ~ 3A^p ~ (-15 ± 6)%. As for the branching ratio, we obtain BR(B° cpK*'^) ~ 
(1.33±0.25) X 10^^, in good agreement with the world average (9.5±0.9) x 10~^ |ly], while without 
NP corrections the result becomes a much smaller value of ~ 5.8 x 10^^. 

B. The second scenario w^ith Oi7,i8,25,26 absent 

In the second scenario, the NP is characterized by O15, 16,23,24 operators and the only relevant 
NP parameter is 023 = \a23\e'^'f'^^ e^^^^ with 4>23 and 823 being the NP weak and strong phases, 
respectively. Following the same way as in the first scenario, we show the results in Tables and 
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TABLE VI: Comparison between the NP predictions and data The NP related observables 

are denoted by "(*)". The second errors for BRs come from the uncertainties of input parameters, 
and the first ones are obtained with the constraint < y?,:„ + 1. The world average for total BR 



BABAR 



BELLE 



The 1st scenario 



The 2nd scenario 



-Soo/Aoo 
-^\\\\/^\\\\ 

S||||+Sxx 



-0.22 ± 0.10 
0.04 ± 0.08 



0.50 ±0.07 
0.25 ±0.07 
0.25 ±0.07 
-0.39 ±0.14 
-0.49 ±0.14 
-0.09 ±0.10 
0.07 ±0.12 
0.02 ±0.10 
-0.09 ±0.06 
0.10 ±0.06 
-0.01 ±0.06 
-0.11 ±0.14 



0.53 ±0.04 
0.21 ±0.02 
0.26 ±0.02 
-0.33 ±0.04 
-0.44 ±0.05 
-0.19 ±0.01 
-0.05 ±0.04 
-0 01+°"" 
0.00 ±0.00 



oi+o-oo 
-0.01 ±0.01 
0.06 ±0.06 



(9.0+°;^ ± 1.9) X io-« 
(4.7 ±0.1 ± 1.0) X 10-^ 
(1.9 ±0.3 ±0.4) X lO"** 
(2.4 ±0.3 ±0.5) X 10-'' 

-0.01 ±0.01 

-0.02 ±0.02 

-0.05 ± 0.05 
0.03 ± 0.03 

-0.01 ±0.01 



0.51 ±0.04 
0.26 ±0.02 
0.23 ±0.02 
-0.49 ± 0.07 
-0.49 ± 0.07 
-0.01 ±0.00 
-0.05 ±0.06 
-0 00+° "° 
0.00 ± 0.00 
-0.00 ±0.00 
-0.00 ± 0.00 
0.05 ±0.07 



(8.7+°-j ± 1.8) X 10-« 
(4.5 ±0.1 ±0.9) X lO-'^ 
(2.2 ±0.3 ±0.5) X lO-'^ 
(2.0 ±0.3 ±0.5) X lO-'^ 

-0.01 ±0.02 

-0.02 ±0.02 

-0.01 ±0.02 

-0.01 ±0.01 

-0.01 ±0.02 



-0.01 ± 0.09 
-0.06 ± 0.10 

-0.10 ±0.25 



IVII where Xmin/'^-'^-f- 0.56/5. The NP parameters in this scenario are given by 

laasi = {l-70tom) x 10~^ ^^23 = 2.36 ± 0.10, = 0.14 ± 0.09, (60) 

with phases in radians. 0,23 produces sizable contributions to the transverse amplitudes. A|||| -|- 
A_L^ ~ Aqo can be understood by following the analysis given in the first scenario. In this scenario, 
because the two terms in both amplitudes Ay and A± in Eq. contribute constructively, we find 
A||||/A^± ~ (^11 /A^ )2 = 1.1. 

As for 4>23 = 0.14 ± 0.09 [= (8 ± 5)°], we obtain = 0.27 ± 0.04,^4° = 0.22 ± 0.04, and 
accordingly S^q = —0.49 it 0.07, A_|_o = —0.05 it 0.06. Since the numerical analysis gives the 
triple-products j4|,,^|i ~ 0.00 ~ 0.01, we therefore obtain S^y = —A\' — A^j^ ~ —0.01 and A^y = 
A^^ — A^j, ~ 0. Note that A_|_o _|_|| are CP-violating observables. We get Ayg = —0.49 ± 0.07, while 
the SM result is Ayo ~ 0.30. For NP related observables, we obtain Syo = -A_lo = 0.05 ± 0.07 
but A^ii T,xx ~ which are rather small. Larger magnitudes of A_|_o and Smq are implied for 



^ It may be better to rewrite as 0,23 = —\a23\e^'^'^''e'^^^^ , where the redefined strong phase is S23 — S23 — tt = 
—0.78 ± 0.10[= (—45 ± 6)°]. The reason is that it is hard to have a large strong phase in the perturbation 
calculation. 
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a larger \4>23\- The BABAR results, displaying arg{A± — A±) ^ and arg(j4|| — ^ 0, may 
hint at the existence of the NP weak phase; consequently, if taking alone the BABAR data, the 
numerical analysis yields (/)23 = 0.23l[J:}^ such that ^^(= 0.29 ± 0.04) / A^(= 0.14j:[];J°), which 
can be rewritten as I;_lo(= - ^t) = -0.43+|];?^ and A^o(= - ^t) = -0.16l|];J^, and 

S||o(= Re{A\\A*o - A\\AI)) = 0.15 ± 0.09. Finally, we get BR(5° (t)K*^) ~ (1.22 ± 0.24) x 10"^ 
which is in good agreement with the world average (9.5 ±0.9) X 10-*^ JJ,]. 

V. SUMMARY AND CONCLUSION 

The large transverse polarization anomaly in the B (f)K* decays has been observed by BABAR 
and BELLE. We resort to the new physics for seeking the possible resolutions. We have analyzed 
all possible new-physics four-quark operators. Following the analysis for the helicities of quarks 
arising from various four-quark operators in the B decays, we have found that there are two classes 
of operators which could offer resolutions to the B (j)K* polarization anomaly. The first class is 
made of Oi7,i8 and 025,26 operators with structures (1 — 75) (8> (1 — 75) and a{l — 75) o"(l — 75), 

respectively. These operators contribute to different helicity amplitudes as i^oo '■ H : ~ 

0{l/mh) ■ 0{l/ml) : 0(1). The second class consists of Ois^ig and 023,24 operators with structures 
(1 + 75) (8) (1 + 75) and a{l + 75) it(1 + 75), respectively, and the resulting amplitudes are given 
as Hoo : H++ : H— ~ 0{l/mb) : 0{l/ml) : 0(1). Moreover, we have shown in Eq. ()33() that 
by Fierz transformation Oi7^i8 can be rewritten in terms of 025,26; and Ois^ig in terms of 023,24- 
For each class of new physics, we have found that all new physics effects can be lumped into a sole 
parameter: 025 (or 023) in the first (or second) class. Our conclusions are as follows: 

1. Two possible experimental results of polarization phases, arg(^_|_) — arg(^||) ~ vr or 0, orig- 
inating from the phase ambiguity in data, could be separately accounted for by our two 
new-physics scenarios with the presence of a large(r) strong phase, 825 (or ^23), and a small 
weak phase, (j)25 (or </>23). In the fist scenario only the effective coefficient 025 is relevant, 
which is related to Oi7^i8,25,26 operators such that ^ H , while in the second sce- 
nario only the effective coefficient 023 is relevant, which is associated with Oi5^i6,23,24 oper- 
ators such that H ^ Note that if simultaneously considering the six parameters 

1^251) ^25, </'25, |«23|, <^23, 023 in the fit, the final results still converge to the above two scenarios. 

2. We obtain A|||| ~ 0.8A_|__|_ in the first scenario, but A|||| ^ A_|_^ in the second scenario. 

3. Our numerical analysis yields ~ 0.10 and ~ —0.19 in the first scenario, but 
gives ~ 0.01 and S^y ~ —0.01 in the second scenario. These two scenarios can 
thus be distinguished. Furthermore, a larger magnitude of the weak phase, <j)25 or ^23, can 
result in sizable A_|_o,S||o- As displayed in Table IVTl we obtain A_|_o ^ ~^||o — —0.05 for 

</'25,(23) = -0.11 (0.14). 

4. The NP related observations 5]qoj|||^^_l, A_|_|| are only marginally affected by weak phases 

</'25,23- 
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5. We obtain BR(S (p^*) — (1-3 ± 0-3) x 10 ^ in two scenarios. Note that we have used the 



rescaled LCSR form factors in Ref. | 2C 
in explaining B K*^, Xg'j data [2J 



2J], where smaller values for form factors were used 
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APPENDIX A 

The LCDAs of the vector meson relevant for the present study are given by [s^l 



(Al) 

(y(P',e)|gi(y)7^75g2(x)|0) 

= -fy(l- !Llhl±^\mve.,p,e*yPz- ei(«p'-W-) (A2) 

V /y "^v / Jo 4 

{V{P\ e)\qM^,uq2{xm = -if^ C e^(«P'-«'-) {el,p', - e^^p'^) ^^{u), (A3) 



where z = y — x with = 0, and we have introduced the light-like vector p'^ = P'^ — rnyZp^/ {2P' ■ z) 
polarization vectors are defined as 



with the meson's momentum P'^ = niy. Here the longitudinal and transverse projections of the 



= 4m = <-4.- (A4) 

Note that these are not exactly the polarization vectors of the vector meson. In the QCDF calcu- 
lation, the LCDAs of the meson appear in the following way 

ly{P',e)\qi^{y)q25{xm = \ dn e*("f'-«'-^) 

4 JO 

+ . (A5) 
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Note that to perform the calculation in the momentum space, we first represent the above equation 
in terms of z-independent variables, P' and e*. Then, the light-cone projection operator of a light 
vector meson in the momentum space reads 



Ml 



rV 



with the longitudinal projector 



MY = ^^^^^^^-^ri_$„(^: 



(A6) 



(A7) 



k=up' 



and the transverse projector 



+ 



fvmv 



'dv{<!>\\{v)-gi\v)) 



6'e* 



_d 



p 

n_n 



du 



P 



dk_i_^ 



k=up' 



(A8) 



where = (1, 0, 0, — 1), n^|^ = (1,0,0,1), is the transverse momentum of the qi quark in the 
vector meson, and the polarization vectors of the vector meson are 



e • nj 



In the present study, we only consider the leading contribution in AQCT)/nT-b for . In Eqs. ()A1|) . 
(|A2|) and HA3|) . <I>||,<I>^ are twist-2 LCDAs, while gj^\g^f^ are twist-3 ones. Applying the equation 
of motions to LCDAs, one can obtain the following Wandzura-Wilczek relations 



#(n) 



u 



V 



dv + 
dv + u 



1 ^\\(v) 



■ dv 



V 

1 $||(u) 

V 



dv 



+ ... 



(AlO) 
(All) 



10 V 

where the ellipses in Eqs. HAlOf) and I|A11I) denote additional contributions from three-particle 
distribution amplitudes containing gluons and terms proportional to light quark masses, which we 
do not consider here. Eqs. HA10|) and (|A11|) further give 



+ 9±{u) 



du 



^ll{v) 



dv + ... , 



{^\\{v)-gi\v))dv 



" $11 (f) 

V 



dv — u 



<^u{v) 



dv 



+ 



(A12) 
(A13) 



After considering Eqs. I|A10|) . (|A11I) . (|A12|) and HA13|) . in Eq. are actually equal to zero. 
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